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A new family of dynamic states are found in a stack of inductively coupled intrinsic Josephson 
junctions in the absence of an external magnetic field. In this state, {2mi + 1)ti phase kinks with 
integers m;'s stack along the c axis and lock neighboring junction together. Large dc power is 
pumped into plasma oscillation via kinks at the cavity resonance. The plasma oscillation is uniform 
along the c axis with the frequency satisfying the ac Josephson relation. Thus this state supports 
strong terahertz radiation and seems to be compatible with the recent experimental observations. 
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The discovery of the intrinsic Josephson effect [ll in 
highly anisotropic layered high-Tc superconductors has 
opened up a new direction for generation of terahertz 
electromagnetic (EM) waves, which have many promising 
applications in materials science, biology, security check 
and so on[2,3. Much effort has been taken to investigate 
the feasibility of such a technique. One idea is to use 
motion of Josephson vortices lattice to excite Josephson 
plasma, and it attains certain degree of success [J [^, 0, 
[3,3. In 2007, it was demonstrated surprisingly that a 
coherent terahertz EM waves were emitted from a mesa 
of Bi2Sr2CaCu208+5(BSCCO) single crystal even in the 
absence of magnetic field [9,]. The experiments [^, 
raise several questions on the mechanism of radiation. 
First, although the working principle is certainly the ac 
Josephson effect, it is not transparent how dc input power 
is pumped into ac plasma oscillation, rather than being 
dissipated as Joule heating. Secondly, how can in-phase 
plasma oscillation can be attained in all junctions? 

The experiments attracted many theoretical 



attentions [ll|, [ij, [13|, [lj|. A new dynamic state is 
proposed to explain the experiments [ill, [l^ . In this 
state, there are (2m/ + l)7r kinks localized at the center 
of the mesa and stacked periodically along the c axis. 
The kinks help to pump large energy into plasma 
oscillation, and a part of which is emitted outside. The 
kink state has already captured the key observations 
of experiments, and is believed to be relevant to the 
emission in experiments. 

The inductively coupled sine-Gordon equation appro- 
priately describes the dynamics of gauge invariant phase 
difference P/ in a stack of intrinsic Josephson junc- 
tions (I J Js) 

dlPi = (1 - CA(2))(sinP/ + pdtPi + P/ - Jext), (1) 

where Pi is the gauge- invariant phase difference at the /th 
junction, A(2)Q/ = Q/+i + Q/_i - 2Q/, ( = A2^/(s + D f 
the inductive coupling, /? the normalized conductance. 



Joxt the bias current. Aab is the penetration depth along 
the c axis and s{D) is the thickness of superconduct- 
ing (insulating) layer. Length is normalized by penetra- 
tion depth in the ah plane Ac and time is normalized 



by Josephson plasma frequency [ll|. The typical value of 
(i is 0.02 and C is 7.1 x 10^ for BSCCO. Here we have 
neglected the variation of the amplitude of order param- 
eter. The treatment is reasonable because the Josephson 
coupling is extremely small in comparison to the conden- 
sation energy. The change in P; , which is associated with 
that in Josephson energy, should not cause appreciable 
change in the amplitude, which is associated with the 
condensation energy. 

In experiments, the thickness of a stack of IJ Js is much 
smaller than the wavelength of electromagnetic wave, 
which makes the EM wave transmission from IJJs to out- 
side very difficult. Same situation happens in a thin ca- 
pacitor. Such a phenomenon is knows as impedance mis- 
match. The mismatch in impedance was first formulated 
explicitly by L. N. Bulaeviskii and A. E. Koshelev[l5[, 
and then confirmed by experiments [l^ with a cylindrical 
As a good approximation, we may neglect 
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mesa 

the effect of radiation on the dynamics of phase. In the 
present paper, we will use the non-radiating boundary 
condition dxPi = 0. A self-consistent treatment will be 
reported elsewhere. 

Equation ([l]) allows for a variety of dynamic states, 
such as state without kink, state with kink and state 
with solitons, each of which has unique electrodynamics 
properties [l9|. The state with kink[ll| is described by 

Pi[x,t) ^ uot + Pi" [x) +Re[-ig{x) e-x]i{iujt)], (2) 

where the firs term at the r.h.s. is the rotating phase with 
frequency uj, the second term the static kink and the last 
term the plasma oscillation whose frequency obeys the 
ac Josephson relation. The plasma oscillation is uniform 
in all junctions which can be realized in a thick stack of 
IJJs. For the first cavity mode, g{x) = Aicos/cix with 
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FIG. 1: (Color online). Several examples of periodic arrange- 
ment of P" in a stack of IJJs calculated from Eq. Here 
only even period is shown. 



FIG. 2: (Color online). Several examples of periodic arrange- 
ment of Pi" in a stack of IJJs calculated from Eq. Here 
only odd period is shown. 



kj = jn/L^. It will be shown later that g(x) contains 
other modes, such as feg, A:2, ^3... as well, but it is sufficient 
to take only fci mode when Ai < 1. The solution with 
kink in Eq. ^ is not written artificially. Actually, it is 
observed repeatedly by computer simulations which 
strongly suggests that this solution is very stable. 

For ease of theoretically treatment, we consider the 
region with < 1. Substituting Eq. ^ into Eq. (P) 
and linearizing the plasma part[l9|, we have the equation 
for the static kink Pf 



cos(fcix)A(^) exp(-iPf). 



(3) 



Equation ([3]) has (2m/ + l)7r kink solutions. The operator 
A^^^ allows for a variety of arrangements of Pj^ along the 
c axis. 

Here we consider several typical periodic arrange- 



ments, i.e. P? 
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with T the period. The results 



for T = 2,4,6 are shown in Fig. [T] (2mi -|- l)7r kinks 
are stacked periodically along the c axis, which takes the 
advantage of huge inductive coupling to lock junction to- 
gether. The kink runs sharply from to (2TOi-|-l)7r in the 
region of width Ap = l/-\/Cl^i|- tti6 BSCCO system, 
C « 10^ which renders the kink almost a step function. 
The huge inductive coupling also makes the kink very 
rigid. As shown in Fig. [Tl even for the same period, 
there are many different ways to pile up the kinks. Thus 
the kink state occupies finite volume in the phase space, 
which makes the state easy to access. The phase kink 
for higher cavity modes can be straightforwardly con- 
structed from the fundamental one. The kink is always 



at the nodes of oscillating electric field. 

With the relation between magnetic field and phase 
(1 — CA(^^)i3/ = dxPi, it is easy to show that flux asso- 
ciated with the static phase kink Pf is (1 — (A'-^^)^i — 
(1 - CA(2))D J^- B^dx = D{2mi + l)7r. Thus the static 
vortices is quantized, depending on Pf and the arrange- 
ment of P( . 

Since no external magnetic field is applied, the total dc 
flux is zero in the IJJs, i.e. = 0. Thus the value of 

kink summed over all junctions vanishes, i.e. + 
1) — 0. For even period such as t = 2, 4, 6 shown in Fig. 
[2 the (2m; -I- l)7r kinks commensurate with the period 
very well. But for odd period, there is no way to put 
(2m/ -|- l)7r kinks in junctions while keeping total dc flux 
zero, which makes the kinks incommensurate with IJJs. 
As depicted in Fig. [2l the kinks adjust themselves in 
order to eliminate the static flux. The resulting kinks 
deviate slightly from (2m; -I- 1)tt. 

From the equation for the frequency component uit, we 
obtain the expression for Ai 



with 



Ai 



-2i 



ikf — iixP' — Plu 



(4) 



Fi = — — / exp(iP;^) cos(fcix)da; (5) 
Lx J 


Fi represents the coupling between the cavity mode and 
plasma. Here only the mode ki of exp(iP*) is taken. In 
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FIG. 3: (Color online), (a) Amplitude of plasma oscillation 
|Ap and (b) IV characteristics. 



principle, there are many different modes fcj, and we have 
to write it down in g{x) at the very beginning. However 
only the mode ki is dominant when the plasma oscillation 
is not strong. The static kink Pf plays a role of pumping 
large dc energy into plasma oscillation. This becomes 
more transparent if we look at the total phase. The total 
phase is rotating with frequency lo but it has kink over the 
X direction. This type of spatial structure is compatible 
with the oscillating electric field, that is the center of the 
kink coincides with the nodes of electric field. As the 
phase is rotating for the time being, the electric field is 
excited. 



Adopting the step- function approximation for Pf 17|, 
Ai can be readily evaluated at different uj. The ampli- 
tude at other cavity modes Aj can be easily generalized 
from Eq. ([4]). The results are also depicted in Fig. [3l 
There is a big enhancement in the amplitude at each 
cavity resonance, which makes the higher frequency har- 
monics visible in the frequency spectrum as calculated by 
computer simulation The higher frequency harmon- 
ics correspond one by one to the high modes in exp(iP;''). 

The radiation power can be calculated by using an 
effective impedance Z (for simplicity, Z is taken as 
real). When Z is large, the radiation can be treated 
as perturbation, and the radiation power is give by 



S = 



\Ai\ /2Z. The dependence of S* on Z is non- 



trivial. For a large Z, much more power can be pumped 
into plasma oscillation, and \Ai\ is larger than that of 
smaller Z, which makes the radiation power large. In 
Ref. [llj, by solving Eq. ([T]) with radiating boundary 
condition numerically, the maximum power at the first 
cavity mode is about lOmW with a mesa of size simi- 
lar to the experiments [§| even for Z = 1000. The power 
estimated by theoretical calculations is still much larger 
than that in experiments. There is an optimal value Z 
which presumes the maximum radiation power (20| . 

For theoretical tractability, we have neglected the ef- 
fect of radiation on the phase dynamics. It does not mean 
that the radiation will kill the kink state. Actually, it was 
found that the kink state survives even for Z = 10, but 



the height of current steps is reduced because it is diffi- 
cult to excite cavity resonance with small cavity quality 
factor 11 1 . 

The IV characteristics is given by the current conser- 
vation law 



Joxt = (3uj + {sin Pi) = fiuj 



/3w|Fi|V4 



(fc2 - t^2)2 + ^2^2 



(6) 



where (• ■ ■)xt represents the average over space and time. 
As considerable amount of dc power is converted into 
plasma oscillation at cavity resonances, the IV deviates 
significantly from linear ohmic behavior and self-induced 
current steps appear, as in Fig. [3] In the absence of ir- 
radiation, the input power Jcxf^ is consumed by dissipa- 
tion by normal current Pui"^ and by dissipation in plasma 
oscillation /3w|Aip/4. 

In conclusion, {2mi + l)7r kink state is formulated in 
the present paper. In this state, the plasma oscillation 
is uniform along the c axis, while the static phase kinks 
of {2mi + 1)tt sit periodically along the c axis. The kink 
induces a new type of current steps in the IV character- 
istics. This state supports intensive terahertz radiations. 
The recent experiments on THz radiations from BSCCO 
single crystals can be interpreted in terms of this state in 
a consistent way. 
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